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ABSTRACT 

Effects of aircraft dynamic characteristics on passenger ride 
quality are investigated to determine ride-quality isocontours similar 
to aircraft handling-qualities contours. Measurements are made on a 
moving-base simulator while varying the aircraft short-period arid 
Dutch Roll frequencies and dampings. Both pilot ratings and subjective 
ride-quality ratings are obtained during "flight." Ride and handling 
qualities were found to be complimentary for the Dutch Roll mode, but 
not for the short-period mode. Regions of optimal ride and handling 
qualities are defined for the short-period mode, and the effects of 


turbulence levels studied. 
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CHAPTER I 








INTRODUCTION AND PROGRAM DESCRIPTION 


Up until the present time, there have been two completely separate 
and independent efforts to analyze and improve the ride quality and 
handling qualities of aircraft. The subject of this dissertation is to 
define, in a quantitative way, the relationships between these two 
quantities and to determine what trade-offs exist between the two. 

This is accomplished by varying the dynamic characteristics of a 
simulated aircraft in a motion-based simulator and obtaining simulta- 
neous measurements of ride qualities from test passengers riding aboard 
the simulator and handling qualities from a pilot flying the simulator, 
Predicted levels of ride quality are also computed using comfort models 
developed by the University of Virginia based on motion parameters of 
the simulator and the particular aircraft design being simulated. 
Finally, isocontours for both handling qualities and ride qualities 
are defined enabling future designers of transport aircraft to weigh 
the merits of their designs from both the passenger's and pilot's 
viewpoints. In addition, this can help a designer weigh the relative 
effects which some modification to an existing design would have on 
ride quality and passenger acceptance as well as handling quality and 
pilot acceptance. Also, the limitations of using ground-based 
simulators with restricted motion capabilities to measure ride qua] ity 
are demonstrated. 

For many years, a considerable amount of work, resources, and 
time have been spent by military and civil groups in the aeronautical 


field to determine what qualities or characteristics an aircraft should 
have for it to be most easily and effectively flown. As a result, 
certain basic criteria have been formulated as guidelines for satis- 
factory aircraft handling qualities (l)-(8). In general, these 
criteria have been determined using variable stability aircraft. 
Stability parameters are varied within a limited range and the opinions 
of numerous pilots on the handling qual ities of each configuration 
obtained. Handling qualities measurements are typically made in 
reasonably calm air, and results thus obtained are generally accepted 
to apply to high turbulence cases as well. This is not necessarily a 
prudent course to follow, in that turbulence can indeed have an effect 
on pj lot opinion. 

Gradually, the physical quantities which have a bearing on 
handling qualities were defined and the range of these parameters 
which delineate satisfactory and unsatisfactory handling qualities 
determined. Naturally, aircraft with different missions have vastly 
different requirements , thus handling qualities are different for a 
bomber or transport- type aircraft than for a high performance, highly 
maneuverable fighter. 

However, even though much attention has been given in the past 
to developing aircraft with good handl ing qual i ties , very little 
attention has been devoted to developing aircraft with good ride 
qualities, i.e., aircraft in which the traveling public find flying 
pleasant (9). In the recent past, with the evolution of high-altitude 
jet aircraft, this has not presented serious problems since passengers 
on commercial flights, as a whole, were not especially annoyed by the 


experience of flying. However, with the increasing use of short 
take-off and landing aircraft (STOL) in conmercial operation, the 
question of acceptable levels of ride quality has arisen due to the 
often unpleasant nature of motion encountered on such aircraft (10). 

STOL aircraft often exhibit objectionable ride qualities due to 
lower wing loadings or higher lift coefficients than aircraft used 
in conventional operations. Designs which use propulsive lift or 
thrust vectoring to attain STOL operations are generally not subject 
to ride-quality degradation because during the cruise portion of 
flight, these craft operate as high-wing loading conventional aircraft 
It is not these designs at which this study is directed. Rather, it 
is the designs which achieve STOL performance by using low-wing 
loading. Wing loading is a measure ■ot the lifting ability of a wing 
per unit of surface area. (For equilibrium flight, it is W/S, the 
aircraft weight divided by the wing reference area.) The lower the 
wing loading, the more susceptible the configuration is to external 
disturbances. Thus, low-wing loading aircraft are more severely 
disturbed by flight through rough air than high-wing loading designs, 
so it is to be expected that STOL aircraft would have worse rids 
qualities than conventional aircraft. It is just such aircraft with 
low-wing loadings which require modifications to improve their ride 
qual ities. 

During the past decade, considerable research has been conducted 
to determine what environmental and psychological factors define a 
person's state of comfort or discomfort. More specifically, human 



response to motions of various forms, noise, temperature, and pressure 
levels have been measured and tolerance levels have been defined (11). A 
recent survey of past work done in this area is contained in Reference 
(12). With regard to the environment experienced by the flying public, 
the University of Virginia has been involved since 1970 in defining 
what factors are involved in determining passenger comfort, and have 
developed several comfort models based on aircraft motion param- 
eters (13). This University of Virginia effort is based on 

simultaneous in-flight measurements of aircraft motions and sampling 
of passenger opinion on regional and commuter airlines on the East 
Coast. 

To improve the ride quality of STOL aircraft, several means have 
been investigated (14). In general, these methods consist of placing 
sensors in the aircraft which sense aircraft motion, usually linear 
accelerations and angular rates. These signals are then used to 
deflect control surfaces which generate aerodynamic forces and moments 
which tend to minimize the motion which the passenger feels. One of 
the disadvantages of some of these systems is that they may tend to 
degrade the handling qualities or controllability of the airplane, 
making it more difficult or annoying for the pilot to fly. In addition 
to a weight penalty such systems may impose on a design, the failure 
of such a ride-control system may present a safety hazard by severely 
increasing pilot work load due to the corresponding change in handling 
qualities, or by exceeding design structural limits. 

Rather than using active control systems to control ride quality, 
one might possibly des Ign a ircraft so that they are inherently pleasant 


to ride. Thus, the purpose of this study is to determine the relation- 
ship between characteristic aircraft motions and aircraft ride quality. 

Most aircraft have five distinct characteristic motions, two 
longitudinal and three lateral. These motions are determined by 
aircraft geometry, mass distribution, and flight conditions such as 
velocity and air density. The phugoid longitudinal mode and spiral 
lateral mode are normally of such long period that these pure motions 
would normally not be sensed by flying passengers. Periods and times 
to double/half amplitude of 30 seconds to two minutes are common for 
these modes. !n fact, these modes are rarely seen in typical flight 
because these motions are readily damped out (usually unconsciously) 
by the pilot. Likewise, the rolling mode is not deemed important to 
aircraft ride quality because of the pilot's tendency to keep the wings 
level in cruise, and when maneuvering, to keep rolling rates small 
(usually less than 10 deg/sec). The two remaining aircraft modes, the 
Dutch Roll and the short-period modes, are of particular interest in 
ride-quality studies since their associated periods and amplitudes 
fall into the spectrum of motions found most uncomfortable by human 
beings (0 - 20 Hz) (12). The quantities which usually define the 
handling qualities of these two modes are the Undamped natural 

frequency, go , and the damping ratio, C , of the short-period mode, 
n s 5 

and the number of cycles to half amplitude, C, , time to half amplitude, 

i *d 

T. , and a rol l-to-sidesl ip parameter, [<£/v j ,, for the Dutch Roll mode. 
* 5 d e d 

Using the parameters established for defining satisfactory handling 
qualities for these two aircraft motions, the limits which satisfactory 
ride quality place on these parameters will be determined by subjecting 


human subjects to such motion in aircraft simulators and eliciting 
their subjective comfort responses, and by using comfort models based 
on computed motion parameters for the simulator and the aircraft 
designs being studied. 

The test program is divided into two distinct phases. The first 
phase investigated feasibility and the effects of varying certain 
parameters on ride and handling qualities. The range of parameter 
variation and the effects of these variations on ride quality were 
studied in the University of Virginia's Analog Flight Simulator. 

Once these studies were completed, the second phase was initiated 
at NASA's Langley Research Center. Here tests were conducted on the 
Visual Motion Simulator (VMS) using aircraft parameters determined to 
be important in the first phase. Simultaneous measurements of both 
ride and handling qualities were made for various aircraft configura- 
tions and finally, the trade-offs between ride and handling qualities 
defined. 
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CHAPTER II 


SIMULATOR EXPERIMENTS AT THE UNIVERSITY OF VIRGINIA 


The University of Virginia's fixed-base analog flight simulator 
(Figure 1 and Appendix A) was programmed with the s ix-degree-of -freedom 
equations of motion (Appendix B) given in Figure 2. The aircraft used 
in the simulation was a 51152 N (11,500 pound) Canadian deHavilland 
DHC-6 Twin Otter, This particular aircraft was chosen because it is 
a typical STOL aircraft and has been in service since 1966 in many 
roles. Also, its flying characteristics are well known, and there 
are many pilots available with flying experience in the Twin Otter to 
validate the ground-based simulations. Flight conditions of level 
flight at 914.4 m (3000 feet) and an equilibrium flight speed of 
78,2 m/s (175 mph) were chosen as the typical environment in which 
this aircraft is operated in present short-haul commuter service on 

the East Coast, The wing loading of the aircraft in this configuration 

2 2 

is 405.2 N/m (27.4 lb/ft ). Based on these flight conditions, 
stability derivatives were obtained from an unpubl i shed NASA document 
containing a mathematical model for the Twin Otter used in a fixed-base 
simulation at the Langley Research Center to study STOL air traffic 
control procedures. These stability derivatives agree well with ones 
contained in Reference (14) for a Twin Otter in approximately the same 
flight conditions. The stability derivatives for this flight condition 
may be found in Table I using standard NACA notation (see Appendices C 
and 0). This condition and its corresponding set of stability deriv- 
atives will be referred to as the "nominal" conditions. 
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table I 


FLIGHT CONDITIONS (EQUILIBRIUM) 


h = 914.4 m (3000 ft) (level flight) I = 55031 kg-m' 


W 


W T 

h* 


- 51152 N (11500 lb) 

= 78.2 m/s (256.67 ft/sec, 175 mph) 
« 1.122 kg/m 3 (0.002177 slug/ft 3 ) 

- 0,045 

= 0.2 ! 


I = 22907 kg-m 2 (16900 slug-ft 2 ) 
I - 37411 kg-m 2 (27600 slug-ft 2 ) 


I =* 1898 kg-m 


xz 


a ( 

c 

b 

5 

6 


(40600 slug-ft ) 

2 

1 o 

(1400 slug-ft ) 
-1.3° 

1.98 m (6.5 ft) 
19.8 m (65 ft) 
39.0 m 2 (420 ft 2 ) 


flap 


= 0 


TURBULENCE CONDITIONS 


a - c = 0.914 m/s (3 fps) 
v w 


a =0 

u 


'M 


a 


a 


m 

a 

* 

'L. 

a 


NOMINAL STABILITY DERIVATIVES 
Long 1 tud 1 na 1 Per i vat 1 yes 

0.3818 C D - 0 

a 

0.045 C = -5.9 

m. 

a 

0.035 C, = 5.504 

q 

5,7295 C D - 0 

' q 

0.1432 C m = -23.948 

-1.9098 C x = -2 C D 

1.52 C - -2 C. 


C - C. - C 
x L 
a 

C = -C. - 

z. L 

a a 

C * -C h 

a : a 

C != -C, 
z. L. 

a a 

C = -C n 
x„ 

q q 

c = -c. 

L_ 

q q 
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TABLE I (Continued) 



-0.89 

- 0.12 

0.1215 


Lateral Derivatives 
C =-0.1 

V 

C. - -0.5488 

C - 0.006 
n_ 



0.5 

0.13 

-0.1855 


Control Derivatives 



0.39 

0.00348 

0.2055 

0.0398 



- 0.1 

- 0.01 

- 1.79 

0.45 
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In addition to the nominal configuration, other stability 
conditions were run. These conditions were produced by varying 

a 

(= 3C /3a. the slope of the pitching moment curve) and C (= 8C /5q * 

' m . m q m 

3C m /36, the pitch damping derivative) in the longitudinal mode, and 

C (= 3C /3$, the side-force damping derivative) and C (= 3C /33, 

y 3 / y n 3 n 

the static directional derivative or weathercock stability derivative) 

in the lateral mode. The longitudinal parameters were varied holding 

the lateral derivatives at the nominal conditions and the lateral 

derivatives were varied holding the longitudinal derivatives at their 

nominal values. 

As an approximation of the short-period mode (15), it can be shown 
that the natural frequency, w n , is given by 



~~ W K ' V 

a q 0 a 



and the damping ratio, ? s , by 


2C s“n 


cm 

2S^ 


(-C -c ) 
ni m . 
q a 





where 

c = longitudinal reference length, mean aerodynamic chord 
I = moment of inertia about the pitch axis 
m = aircraft mass 

2 

= dynamic pressure, JjpUg 
Ug - equilibrium flight speed. 


Thus, by judicious choice of C and C a designer could obtain 
* ~ iff nv 

a q 

any set of a) and 5 for the short-period mode he wouid desire, 
n s s 

holding all other parameters constant. The primary factor which 

determines C is the static margin, the distance which the aircraft 
m a 

center of gravity is in front of the center of pressure. For example, 

the designer, by proper selection of center of gravity location, may 

select a desired value of C . The value of the stability derivative 

m a 

C is determined by the size and location of the horizontal tail. 
m q 

Therefore, by choosing a horizontal tail of appropriate size and 

locating it a specified distance behind the center of gravity, the 

designer can obtain a desired value of C . Hence the designer has 

m q 

control over the defining parameters for the short-period mode, but he 

must be aware of how his variations of C.G, location and tail size and 

location might affect other stability derivatives (e.g., C m ) and what 

a 

ramifications such variations would have on cost, complexity, and 
mission performance of a particular configuration. 

Likewise, C and C are dominant factors which determine the 

y e n 3 

damping and frequency of the Dutch Roll mode and are therefore 
Important to the handling qualities of the Dutch Roll mode. For 
simplicity in this study, only variations in C and G were used 

y $ n 3 

to modify Dutch Roll characteristics. However, other stability 

derivatives, such as C 0 and C , have direct effects of Dutch Roll 

6 n r 

characteristics and, at the designers option, he may wish to vary 

these in addition to, or instead of C and C . The values of 

v 3 n 3 

C and C are largely determined by the size location of the 

y e n e 

vertical tall; thus the designer, by giving appropriate attention to 


the matter, may choose a design which will give him values for C , 

y 3 

C , C , and C 0 which will result in good handling qualities in the 

n e n r *(3 

Dutch Roll mode. 

For each set of stability derivatives to be investigated, the 

fourth-degree longitudinal and lateral equations of motion were solved 

for their characteristic roots. From these characteristic values, the 

short-period undamped natural frequency, U) and the damping ratio, C 

n s s 

were computed from the longitudinal roots. Also, the number of cycles 

to half amplitude, C, , the time to half amplitude, T. , and the roll- 
ed 2 d 

to-s ides lip parameter, U/v f° r the Dutch Roll mode were found from 

the lateral roots. (V and T, are related to the undamped natural 

% ^d 

frequency, 0) , and damping ratio, of the Dutch Roll mode by the 

n d d 


following expressions: 
-lnft) 


27T 




i 4>/v | . is related in a complex way to 0) and C, and is found by 
e d *”*d ® 

the relation 


I <f>/v 


e'd U 


= U/e| d 


where 1 4*/ 3 J ^ is found by a modal analysis of the lateral equations 


solved at the Dutch Roll w. and C.. 

n d d 


It Is these quantities which are used to define a design's short- 
period and Dutch Roll handling qualities (l)-(8). For the longitudinal 
short-period mode, the natural frequency and damping ratio describe the 
response of an aircraft to an abrupt pitch change from equilibrium 
flight and its subsequent return to equilibrium conditions. Similarly, 
for the lateral Dutch Roll mode, C, and T describe the motion of an 
aircraft after being disturbed from equilibrium heading, roll, and 
sideslip conditions. The rol 1 -to-s idesl ip parameter is a measure of 
the roll induced for a unit lateral gust and is a gauge of an 
aircraft's lateral responsiveness. 

The analog flight simulator was programmed with various combina- 
tions of stability derivatives to cover, as well as possible, the 
regions in which handling qualities are most often defined for the 
short-period and Dutch Roll aircraft modes. Coupled to the analog 
computer was an electronic noise generator adjusted to disturb the 
simulated aircraft with .914 m/s (3 ft/sec) rms turbulent gusts in the 
normal and lateral directions. This was accomp 1 i s hed by superimposing 
the random electronic signal on the ct and g variables in the analog 
equations of motion. The simulated turbulence was essentially random 
noise with the spectrum shown in Figure 3. 



FIGURE 3. GUST INPUT POWER SPECTRUM 


For each different aircraft configuration, the simulator was 
operated at least 12 times while a pilot flew the simulator attempting 
to maintain straight and level flight for over 200 seconds. This run 
duration was chosen to permit the rms normal and transverse accelera- 
tions computed by the analog computer to stabilize at a steady-state 
value. The associated comfort rating for each flight was found by 
using an empirically-derived comfort model developed at the University 
of Virginia (16) using their comfort- rating scale. (See Table H.) 

This particular model was chosen for its relative simplicity and because 
its defining parameters could be obtained from the analog computer with 
relative ease. The model was derived by simultaneously recording air- 
craft motion and sampling passenger opinion of ride quality on regional 
and commuter airlines during actual flight operations on the East Coast. 
This data was then statistically analyzed and a best fit curve deter- 
mined (13). Finally, the average rating for each aircraft configuration 
was found and converted into a passenger satisfaction level by a 
statistically-determined transformation (17), shown in Figure 4. This 
relationship was formulated by analysis of questionnaire data recorded 
on the above-mentioned commercial flights asking the passengers to 
report their comfort levels and their willingness to take another fl ight 
based on their recent flight environment. 

The values of passenger satisfaction due to variations in the 
short-period handl ing qual i ties are plotted in Figure 5. The solid 
lines in this figure indicate the presently accepted boundaries for 
short-period handling qualities. The dashed lines indicate lines of 
constant ride quality, as suggested by the data points. The trend is 
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TABLE II 

COMFORT RATING SCALE 

1 - Very comfortable 

2 - Comfortable 

3 - Neutral 

4 - Uncomfortable 

5 - Very uncomfortable 


COMFORT MODEL 

CR = 2 + 13.8 a N + 4.52 a T - 2.816 

where a^ - normal rms acceleration (g's) 

a T = transverse rms acceleration (g's) 
CR - comfort rating (1-5) 
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MATHEMA TICAL FIT TO THE CURVE 
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FIGURE 4. RELATIONSHIP BETWEEN PASSENGER SATISFACTION 
AND COMFORT RATING 
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FIGURE 5. SHORT-PERIOD HANDLING QUALITIES 
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FIGURE 6. DUTCH ROLL HANDLING QUALITIES 
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for increasing passenger satisfaction as the damping ratio and the 
undamped natural frequency increase. These trends are to be expected, 
since an increase in damping ratio will directly result in decreasing 
normal acceleration and improving ride quality. Further, increasing 
the natural frequency tends to make a design more responsive, allowing 
a pilot to better maintain control and keep motion to a minimum, also 
improving ride quality. 

The effects of variations of Dutch Roll parameters on ride quality 

• 1 . ' ; ' 

are shown in Figures 6 and 7. The solid line indicates boundaries for 

regions of acceptable Dutch Roll handling. Ride quality and passenger 

acceptance generally improve as , Tj_ , and |<J>/v 1 , decrease, although 

2 d 2 d e 

the trends are not as clear as those of the short-period mode. Again, 
this might be expected since decreasing the number of cycles to half 
amplitude and the time to half amplitude corresponds to increasing 
the damping of the Dutch Roll mode and reducing lateral acceleration. 
Likewise, decreasing the rol 1-to-sidesl ip parameter implies reduction 
of the aircraft roll response to lateral gusts resulting in lower lateral 
accelerations and better ride qualities. Also, it appears that the 
changes in longitudinal short-period parameters had a greater effect on 
ride quality than did the changes in the Dutch Roll parameters. 

No attempt was made in this part of the study to measure pi lot 
opinion since a pilot experienced in evaluating handling qualities was 
not available and no motion or visual cues are provided by the simulator. 


CHAPTER ill 


SIMULATOR EXPERIMENTS AT THE NASA LANGLEY RESEARCH CENTER 

Using the data and experience gained during the tests on the 

University of Virginia's Analog Flight Simulator, 27 different sets 

of stability derivatives were selected to be studied further at the 

NASA Langley Research Center. Fourteen longitudinal cases and 12 

lateral cases were chosen in addition to the "nominal" set of Twin 

Otter stability derivatives contained in Table I of the previous 

section. These cases and the corresponding flight conditions are 

contained in Table III. These cases were chosen to cover the regions 

in which longitudinal and lateral handling qualities are most often 

defined for the short-period and Dutch Roll modes by aircraft 

designers. Figure 8 shows the contours which define longitudinal 

short-period handling qualities as a function of undamped natural 

frequency, to , and damping ratio, 5 . The numbered points on this 
n s S 

figure indicate the 14 test cases studied, plus the nominal case. 
Similarly, Figures 9 and 10 depict the 12 lateral test cases studied, 
plus the nominal case plotted in terms in which the Dutch Roll 
handling qualities are normally defined, namely the number of cycles 
to half amplitude, the time to half amplitude, and the rol 1 -to-s ides 1 ip 
parameter. 

The Visual Motion Simulator (VMS) at the NASA Langley Research 
Center, a synergistic motion-based simulator with the basic interior 
and instrumentation of a jet transport cockpit (Figures 11 and 12), was 
programmed with the 27 test cases described. A schematic diagram of 




TABLE III 


FLIGHT CONDITIONS (EQUILIBRIUM) 
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h = 914.4 m (3000 ft) (level flight) 

W = 51152 N (11500 1b) 

U Q = 78.2 m/s (175 mph = 256.67 ft/sec) 
p =1.122 kg/m 3 
C T =0.045 
h* = 0.2 

I =22907 kg-m 2 (16900 slug-ft 2 ) 

X ■ 

I = 37411 kg-m 2 (27600 slug-ft 2 ) 


xz 

“0 

c 

b 

S 

5 flap 


= 55031 kg-m z 
(40600 slug-ft ) 

= 1898 kg-m 2 „ 
(1400 slug-fr) 

= -1.3° 

= 1.98 m (6.5 ft) 

= 19.8 m (65 ft) 


39.0 m 2 (420 ft 2 ) 


= 0 


TURBULENCE CONDITIONS 
(DRYDEN MODEL) 


Spectra 


<|> (od) = a 2 -r- 
Y U V U TtU 


0 1 + (L u w/U o r 

, M . . > £j 1 1 * M 


0 \ [1 + (L V UJ/U 0 ) 2 ] 2 


2 L w / 1 + 3(L w‘ 

J, w (u) = a w iu- r — ; — 
U 0\[l + (Lji) 


3(L W “/U 0 )' 


(l .oj/u 0 ) 2 3 2 


a u = 0 y - cr w = 1.524 m/s (5 ft/sec) 


L u = L v = L w = 762 m (2500 ft) 


x = 


Pilot Position with Respect to Aircraft c.q. 
2.7 m (8.8 ft) y = -.49 m (-1,6 ft) z = 0 








TABLE HI (Continued) 
LONGITUDINAL TEST CASES 


Case 

Number 

C 

m 

2L 

C 

_!!a 

Short-period Mode 

OJ 

n s 

l£P. s JL ?s 

2 

-1.9098 

0 

0.577 

0.332 

3 

-4.0 

-20.0 

0.885 

0.377 

4 

-0.64 

-23 . 948 

0.463 

0.782 

5 

-6.76 

-23.948 

1.132 

0.320 

6 

-3.24 

7.0 

0.731 

0.193 

7 

-0.64 

7.0 

0.286 

0.499 

8 

-1.9098 

-70.0 

0.796 

0.868 

9 

-4.84 

-56.0 

1.041 

0.568 

10 

-0. 16 

-101.0 

0.677 

1.347 

11 

-2.56 

-119.0 

0.978 

1 .065 

12 

-5.76 

-119.0 v 

1.231 

0.846 

13 

-1.9098 

16.0 

0.514 

0.150 

19 

-3.24 

16.0 

0.704 

0.109 

20 

-6.76 

7.0 

I.072 

0.132 


TABLE III (Continued) 




LATERAL TEST CASES 



Case 

Number 

Vi C 

n B ,,, 

C 

1/C, 

% 

Dutch Rol 1 Mode 
(sec“ l ) 

I't’/Vglj 

ffa/s ) 

14 

0.01583 

0.62 

2.485 

0.260 

1.699 

15 

0.81 

-0.89 

0.999 

0.984 

0.766 

16 

0.1215 

-1.875 

3.076 

1.184 

0.160 

17 

0 

-1.0 

2.384 

0.171 

2.736 

18 

0.01417 

0.86 

1.193 

0.115 

1.914 

21 

0.1215 

0 

1.978 

0.743 

0.165 

22 

0.01167 

0.82 

0.504 

0.045 

2.098 

23 

0.015 

0.74 

1.864 

0. 186 

1.810 

24 

0 

-0.75 

1.533 

0.112 

2.737 

25 

0.25 

3.0 

0.107 

0.055 

0.194 

26 

0 

-0.5 

0.721 

0.053 

2.738 

27 

0.0075 

-0.25 

2.939 

0.243 

2.231 
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FIGURE 0. SHORT-PERIOD HANDLING QUALITIES BOUNDARIES 
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FIGURE 9. DUTCH ROLL HANDL I NG QUALITIES BOUNDARIES 
AND LATERAL TEST CASES* (1/C^ vs. 1 4>/v_ | ) 


VERSE TIME TO HALF 
•LITUDE l/T 1/2 (sec' 1 ) 



SATISFACTORY 

HANDLING 

QUALITIES 


UNACCEPTABLE 

HANDLING 

QUALITIES 


ROLL-TO- SIDESLIP PARAMETER, | <#> / V e | (deg/m/s) 


URE 10. DUTCH ROLL HANDLING QUALITIES BOUNDARIES AND LATERAL TEST CASES 

(1/Ti vs. |4»/v e | ) 



FIGURE II. THE VISUAL MOTION SIMULATOR (VMS) AT THE NASA LANGLEY RESEARCH CENTER 


the simulator, its control systems, and its data output capabilities 
is presented in Figure 13. A CDC-6600 digital computer, used 
exclusively to operate the real-time simulators at NASA Langley, was 
programmed with the aircraft flight conditions, stability derivatives, 
six-degree-of -freedom differential equations of motion, Dryden gust 
model, a routine for computing rms values for 22 parameters, and a 
simulator washout routine. The program integrates the nonlinear six- 
degree-of -freedom equations of motion 32 times a second to describe 
the motion of the real aircraft as a function of time. These values 
are used by the simulator washout routine to determine the position 
of the simulator's six moveable legs as a function of time, providing 
motion sensations to the occupants aboard the simulator. Due to the 
limited displacement of the simulator legs, and the dynamic character- 
istics of the hydraulic actuators, the simulator is not capable of 
producing the magnitude and duration of displacements, velocities, and 
accelerations of the real aircraft. It is the purpose of the washout 
routine to appropriately scale down the predicted motions of the real 
airplane to values which the simulator can produce without exceeding 
one or more of its design limits. The washout routine also attempts 
to drive the simulator legs back to their neutral position following 
a disturbance from equilibrium in anticipation of a future disturbance. 
This "centering" routine is to allow for maximum displacement during 
some future motion, and the motion due directly to this portion of the 
routine is designed to be of such magnitude that it be subliminal to 
the simulator occupants (18). (Time constants of five seconds or 
greater are characteristic of such motions.) 



THE VISUAL MOTION SIMULATOR 













FIGURE 13. BLOCK DIAGRAM OF MOTION-BASE SIMULATOR APPARATUS 
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A detailed description of the physical dimensions and the 
performance specifications of the VMS may be found in References ( 1 8 ) , 
(19), and (20). For all tests the pilot was to perform an instrument 
flight rules ( I FR ) task ar;d was given no visual or 1 'out -the -window" 
cues. This condition was chosen to provide the pilot with a 
representative workload of cockpit duties and, at the same time, 
somewhat simplify the simulation. A ride-quality subject with past 
experience using the U.Va. -developed 5-point ride-quality rating scale 
and actual in-flight ride-quality evaluation experience rode aboard 
the simulator with the pilot to evaluate the ride quality of each 
configuration. Like the pilot, the ride-quality subject was given 
no visual cues, and the flight instruments in front of him were 
covered. The pilot used in the simulations was a fully instrument- 
rated pilot with over 3000 hours of flight time. 

The 27 test cases were run in random order three times each, with 
each test run lasting 16 minutes. Each 16-minute segment consisted of 
the following subsegments. For the first 10 minutes, the pilot is 
instructed to fly straight and level. During this time, the ride- 
quality subject on board the simulator evaluates the ride quality every 
two minutes. For the next two minutes, the pilot executes a two-minute 
turn in which the aircraft changes heading by 180°, descending 304.8 m 
(1000 feet) for the first minute and climbing 304.8 m (1000 feet) to 
the original altitude In the second minute. The ride-quality subject 
evaluates the ride of this two-minute segment. The pilot is then 
asked to execute a second two-minute turn, identical to the first, and 
return to the original aircraft heading. Again, the ride-quality 


subject evaluates the ride of this segment. In the final two minutes 
of each run, the pilot is instructed to separately pulse the elevator, 
aileron, and rudder to enable him to better evaluate the handling 
qualities of the particular configuration being investigated. Simul- 
taneously, the ride-quality subject evaluates the comfort levels of 
the motions produced by each of the control pulses. For the study at 
hand, only the ride-quality evaluations taken during the straight and 
level portions of each run will be used for analysis. Evaluations 
taken during the turns and control pulses are recorded for possible 
future use. Following th® run, the pilot is asked to complete a 
questionnaire rating his ability to maintain straight and level flight, 
and give his opinion of the overall handling qualities of the case 
being studied using the Cooper-Harper rating scale (21). The Cooper- 
Harper scale is a 10-point rating scale used by test pilots to 
quantify pilot opinion of aircraft handling qualities, with 1 being 
most desirable and 10 being almost unflyable. Also, the ride-quality 
subject is asked to give an overall rating of the ride quality of the 
cruise portion of this configuration. Samples of the rati rig forms used 
by the pilot and subjects are presented as Figures 14 and 15. 

During each run, continuous strip-chart recordings are made 
displaying time histories of various aircraft parameters. These 
include the three linear accelerations and three angular rates of the 
aircraft in the body axes, elevator, aileron, and rudder deflections, 
throttle position, altitude, rate of climb, airspeed, and heading, as 
well as a time channel indicating the times when ride-quality responses 
are taken. A sample output is shown in Figure 16. 
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Date Case Number 

Rate the following using the Cooper-Harper scale: 

Your ability to maintain straight and level flight 

Your opinion of overall handling qualities ___________ 

Do you have any particular comments regarding the handling qualities 
of this configuration? 
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FIGURE 14. EVALUATION PILOT'S HANDLING QUALITIES RATING SHEET 
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Date 


Case Number 


What is your overall opinion of the ride quality of this configuration? 


i 1 

•C K 


‘j. | 

| | 

♦4 



Rate the overall ride quality of this configuration 



What is your overall opinion of the ride quality of this configuration? 



Rate the overall ride quality of this configuration 


FI GURE 15. RIDE -QUAL 1 TY TEST SUB JECT 1 S RAT I NG SHEET 
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FIGURE 16. TYPICAL STRIP-CHART OUTPUT OF FLIGHT PARAMETERS 
(LEVEL FLIGHT, a = 1.524 m/sec (5 fps)) 




In addition to these measurements, the main program which controls 
the simulator computes rms values for the three linear accelerations, 
three angular accelerations, three angular rates, and all control 
deflections for each two-minute segment of each run. A sample output 
is shown in Figure 17. Rms values of the motion parameters are 
evaluated for both the values predicted for the real aircraft by the 
six-degree-of-freedom equations of motion, and the values computed to 
drive the simulator after accounting for the washout system. Using 
real aircraft and simulator rms acceleration variables, comfort ratings 
are computed for each two-minute segment by empirically-determined 
l inear comfort models contained in Table IV (22) (23). The particular 
comfort models used in this portion of the study were chosen because 
they were more advanced than the one used in the preliminary study. 

The first model incorporates motion in all six degrees-of-f reedom and 
the second appl ies for cases where lateral motion dominates. Both 
models use rms accelerations as the defining parameters, and use of 
more complex models would not be warranted for this study. Finally, 
an inertial package is placed aboard the simulator to sense and record 
the motion of the simulator itself for comparison with the values of 
the computed parameters and to evaluate the errors between the driving 
signals after washout and the actual motion sensed in the cabin of the 
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FIGURE 17. TYPICAL COMPUTER OUTPUT OF RMS MOTION VARIABLES AND PREDICTED COMFORT RATINGS 










TABLE IV 


COMFORT MODEL 

CR = 1.85 + 11.5 a N + 5.7 a T + 1.0 a L 

+ 0,2 a ROLL + 0,2 a PITCH + 1,5 a YAW (a N > 1.6 a T ) 

CR ■ 1.9 + 8.1 a N + 20.9 a T « (a N <: 1.6 a T ) 

where a.. - normal rms acceleration (g 1 s ) 

a^. = transverse rms acceleration (g 1 s ) 

a^ = longitudinal rms acceleration (g 1 s ) 

a B », , = rms angular acceleration about X body fixed axis 

(radians/sec^) 

a d i rr u = rms angular acceleration about Y body fixed axis 
(rad ians/sec^) 

a u „,, = rms angular acceleration about Z body fixed axis 

YAW (radians/sec 2 ) 

CR = comfort rating where; I - very comfortable 

2 - comfortable 

3 - neutral 

4 - uncomf o r tab 1 e 


5 - very uncomfortable 


CHAPTER IV 


RESULTS AND CONCLUSIONS 

From the data collected on the VMS, several statistical quantities 
were computed for each separate test case, i .e. , each unique set of 
stability derivatives and aircraft handling qualities. Computed were 
means and standard deviations of the comfort ratings of each test case, 
These quantities were found for the comfort ratings predicted for the 
actual aircraft, the ratings predicted based on the actual simulator 
motion, and the comfort ratings elicited from the test subjects riding 
aboard the simulator. These quantities are tabulated in Table V. 

Also, typical power spectra of the three linear accelerations obtained 
from measurements made by the inertial package placed aboard the 
simulator are presented in Figures 1 8 , 19, and 20. The rms acceleration 
quantities obtained from this inertial package agreed within 10% of 
those computed for the simulator after washout by the central computer, 
indicating the computer-generated quantities were valid. Only the 
values computed and recorded during the straight and level portions of 
each simulator run were used when computing these statistics since it is 
the ride quality during the cruise portion of flight which is of primary 
interest in this study. Thus each mean and standard deviation is based 
on 15 data points since each test case was run three different occasions, 
and five r ide-qual i ty measurements were taken during the 10 minutes of 
straight and level flight of each run. 

Comparing the mean comfort ratings for the 27 test cases as 
computed from the equations of motion of the real aircraft with the 


TABLE V 



MEAN AND STANDARD 

DEVIATION 

OF COMFORT RATINGS 

FOR TEST CASES 

(n = 15 ) 



Based on Predicted 

Based 

on 

Actual 


Case 

Real Aircraft 

Motion 

Simulator 

Motion 

Subjective Response 

Number 

CR 

a (CR) 

CR 

cj (CR) 

CR 


1 

2.70 

. 13 

2.03 

.04 

2.00 

.00 

2 

2.97 

.12 

2.13 

.08 

2.07 

.26 

3 

2.71 

.13 

2.11 

.08 

2.00 

.38 

4 

2.91 

.13 

2.10 

.10 

2.00 

.00 

5 

2.53 

.10 

2.08 

.08 

2.27 

.46 

6 

2.88 

.08 

2.07 

.05 

2.13 

.35 

7 

3.60 

.62 

2.06 

.04 

2.00 

.00 

8 

2.64 

.05 

2.09 

.08 

2.00 

.00 

9 

2.60 

.14 

2.11 

.08 

2.00 

.00 

10 

2.80 

.07 

2.03 

.05 

2.07 

.26 

11 

2.63 

.06 

2.09 

.06 

2.00 

.00 

12 

2.57 

.09 

2.07 

.07 

2.13 

.35 

13 

3.61 

.29 

2.16 

.08 

2.40 

.51 

14 

3.04 

.65 

2.34 

.16 

2.33 

.49 

15 

3.39 

.09 

2.17 

.05 

2.87 

.35 

16 

2.93 

.34 

2.10 

.14 

2.27 

.46 

17 

4.03 

.56 

2.57 

.16 

2.33 

.49 

18 

4.15 

.84 

2.78 

.32 

2.40 

.51 

19 

3.09 

.16 

2.05 

.04 

2.40 

.51 

20 

2.67 

.10 

2.05 

.06 

2.73 

.59 

21 

2.90 

.17 

2.10 

.12 

2.40 

.51 

22 

4.34 

.57 

2.91 

.26 

2.33 

.49 

23 

3.44 

.72 

2.48 

.21 

2.20 

.41 

24 

3.97 

.84 

2.63 

.33 

2.27 

.46 

25 

5.00 

.00 

2.28 

.09 

2.47 

.52 

26 

4.23 

.86 

2.72 

.28 

2.47 

.52 

27 

2.79 

.15 

2.22 

.07 

2.27 

.46 
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ratings based on the commands to actuate the simulator motion, one 
finds very little agreement between these two quantities (see Table V) 
These results are due to the physical limitations of the motion-base 
cockpit and the corresponding washout routine which controls the move- 
ment of the simulator's six moveable legs. This implies that the VMS 
was not capable of reproducing the ride quality of a given aircraft 
configuration on a one-to-one basis for the present study. This is, 
of course, due to the washout system and the displacement, velocity, 
and acceleration limits of the simulator. The correlation coefficient 
as computed by 



ktt J. (x r* )( y,-y> 


a a 

x y 


between these two quantities is r = 0.75 however, indicating that 
although the motion is not being exactly reproduced as intended, there 
is a definite relation between the simulator ride quality and that of 
an actual ai rcraft . 

Comparing the mean comfort ratings computed from the commands to 
move the simulator legs with the mean of the comfort ratings recorded 
by the ride-qual ity subjects riding aboard the simulator, one finds 
the agreement for most cases to be quite good (see Table V) . The 
correlation coefficient for this data is r = 0.30. Part of the 
discrepancy between the two sets of mean comfort ratings may be 
accounted for by the fact that the comfort model used to predict 
comfort ratings based on rms accelerations is a continuous function 
while the actual rating scale as defined and used by the test subjects 
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is an integer scale. This also accounts for the low correlation 
coefficient for this data. Also the number of responses used in 
determining each mean was somewhat small (n - 15). In spite of these 
limitations, the agreement between these two suggests that the experi- 
enced ride-quality subjects were indeed responding in the "proper 11 
(predicted) manner to the motions to which they were subjected, and 
further validates the empirical comfort models. This tends to support 
the theory that ground-based simulations for ride-quality studies may 
be of significant value, if the desired motions and environmental 
conditions being simulated are acceptably reproduced in the laboratory. 
Using the previously-described transformation (22) which translates 
a comfort rating to the percent of the population satisfied with the 
ride, the satisfaction levels of the comfort predicted for the real 
aircraft for the cases when the longitudinal short-period handling 
qualities were varied, are plotted in Figure 21, superimposed upon the 
contours of short-period handling qualities, in a similar manner, 
variations in the lateral Dutch Roll, mode handling qualities resulted 
in the satisfaction levels shown in Figures 22 and 23, plotted along 
with contours of accepted handling qualities for the Dutch Roll mode. 

From the data of Figure 21, lines of constant satisfaction level have 
been drawn and are shown in Figure 24 along with the contours of 
constant handling qualities for the longitudinal short-period mode. 

It is recognized from this figure that ride quality is indeed a function 
of the parameters which define short-period handling qualities, namely 

(0 and £ , Thus a trade-off exists between longitudinal handling 
"s s 

qualities and ride quality. Hence, from a design standpoint, the 

49 






INVERSE TIME TO HALF 
AMPLITUDE 1 /T, (sec* 


081 . 8 


085.6 


TURBULENCE 
cr = 1.524 m/s I5fps)rms 


082.3 


SATISFACTORY 

HANDLING 

QUALITIES 


025.0 


78.9 

O 


84.2 


UNACCEPTABLE 

HANDLING 

QUALITIES 


0 - 68.0 

048.4 


43.1 


O 5 ! — 

~" 053.3 
046.2 


ROLL-TO- SIDESLIP PARAMETER, I <#» / V e | (deg/m/s) 


FIGURE 23. PERCENT POPULATION SATISFIED, DUTCH ROLL HANDLING QUALITIES 

(l/T, vs. I d>/v_| ) 


UNACCEPTABLE 

HANDLING 

QUALITIES 



UNACCEPTABLE 

HANDLING 

QUALITIES 


TURBULENCE 
1.524 m/s (5fps)rms 



-.1 1 . 


J 


designer can choose an operating point and know how both the handling 
qualities of the short-period mode and the corresponding level of ride 
quality would compare with some other design point. This would enable 
the designer to weigh the trade-offs between ride and handling qualities 
and help him make an intelligent choice of an operating point for good 
handling qualities and give him an idea of what level of passenger 
satisfaction to expect from his design. 

Similarly, from Figures 22 and 23, one can identify trends of 
improving ride quality as defined by parameters in which lateral Dutch 
Roll handling qualities are defined. While it is difficult to con- 
struct contours of constant ride quality due to the sparsity of data 
available, the previously-observed trends hold. Namely, decreasing 
i ar| d 1 4>/v e | d tends to improve ride quality and raise the 
percentage of the population satisfied with the quality of the ride. 

Thus in the lateral mode, good ride quality is complimentary to good 
handling qualities of the Dutch Roll mode. It is worthy to note here, 
however, that while these regions are complimentary, there is a sig- 
nificant amount of variation in the percentage of passengers satisfied 
over the region where Dutch Roll handling qualities are judged satis- 
factory (68 - 85%). The designer should note here that simply 
designing an aircraft with satisfactory Dutch Roll handling character- 
istics is not all he should strive for from the ride-quality standpoint. 
He should realize from these figures that the further he can move his 
operating point away from the boundary between acceptable and unaccept- 
able handling qualities, the more he can improve passenger satisfaction 
as well as improve handling qualities; he may be able to improve 
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passenger satisfaction by up to 20% of the total passengers carried for 
a design of similar wing loading and equivalent turbulence conditions 
as the test case. It should be re-emphasized that all the measurements 
taken were on one particular design, one flight condition, and one 
turbulence level. One would not expect the values produced by this 
study to be universal, applying to all conditions and turbulence levels; 
however, the general nature of the isocontours for ride quality and the 
other trends outlined should hold regardless of wing loading, turbulence 
level, and flight conditions. 

If one takes the linearized longitudinal equations of motion for 

the free airplane with controls fixed, eliminates the force equation in 

the X body axis direction, and assumes variations in u are negligible, 

one obtains an approximation of the short-period mode (24). From these 

equations, one may develop an expression for the normal acceleration as 

a function of w, the forcing frequency of an imposed sinusoidal gust 

oscillation. Integration over all forcing frequencies, one obtains the 

following relations for the rms normal accelerations as a function of 

the gust magnitude, w^, the short-period undamped natural frequency, 

u , and the short-period damping ratio, r : 
n_ s 




V w q 






where Kj, K^, and k are constants. Derivation of these relations 
may be found in Appendix E. 

These relations, while derived by using the assumptions of controls 
fixed and allowing only sinusoidal variations in gust velocities, 
exhibit the contours found experimentally. Holding £ s constant, these 
equations all predict that as the short-period undamped natural 
frequency is increased, a^ is decreased resulting in an increase in 
passenger satisfaction. Similarly holding constant, a^ is decreased 
as ? s increases below £ =1, predicting an increase in satisfaction 

for £ < 1. Also, a., is increased as £ increases above r = l 
predicting a decrease in satisfaction for £ g > 1 . These trends agree 
with those observed experimentally in Figure 24. 

The effect that variation of the turbulence level could ha^e on 

the contours of constant ride quality may be studied using this 

mathematical model. For a given constant passenger satisfaction 

contour, the corresponding comfort rating may be found using Figure 4. 

The rms acceleration required to produce this level of comfort is 

determined from the simplified comfort model (Table. If) neglecting 

transverse acceleration terms. Since the normal acceleration is 

directly proportional to rms turbulence level, one may find, by an 

inverse method, contours of constant passenger satisfaction for any 

turbulence level. For example, <t = 0.914 and 2.134 m/s (3 and 7 fps) 

result in the following satisfaction levels: 
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% Satisfied 


0.914 m/s 

(3 fps) 

89.6 

88,7 

86.9 

82,0 

1,524 m/s 

(5 fps) 

87.0 

85.0 

80,0 

63.0 

2.134 m/s 

(7 fps) 

83.6 

82.1 

72.0 

44.0 


These results are plotted In Figures 25 and 26. For the contours shown, 
passenger satisfaction varies only 7-6% for a w = -914 m/s (3 fps); how- 
ever, for a = 2.134 m/s (7 fps), passenger satisfaction varies 39-6%. 
w 

Thus, for flight in regions of heavy turbulence, the selection of 
short-period characteristics could markedly affect the degree of 
passenger satisfaction. 

The pilot's evaluations of his ability to maintain straight and 
level flight, and his opinion of overall handling qualities for all test 
cases are tabulated in Table VI. The pilot made several trial runs on 
several test cases before actual testing was begun, and these are also 
included in the table. This was done to let the pilot gain fam iliarity 
with the simulator, its controls and instrumentation, and use of the 
Cooper-Harper pilot-rating scale, The results compiled in the table 
contain the pilot ratings for all test cases listed in the order in 
which they were taken, included are the trial run cases as well as the 
three test cases, during which ride-quality measurements were taken. 
Referring to Table VI, one can see that for most cases, the pilot was 
consistent among the three ratings he used to evaluate a particular 
case, or agreed closely on two of the three evaluations. Only on one 
case, number 25, did the pilot return three substantially different 
ratings for each test run. The final column in Table VI lists the 
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TABLE VI 


PILOT RATINGS (COO PER- HARPER SCALE) 


Case 

Number 

Ab 1 1 ity to Maintain 
Straight and Level 
FI i qht 

Opinion of 
Overal 1 

Handl inq dual i ties 

Predicted 

Cooper-Harper 

Rating 

1 

3.5, 4, 2, 5, 3 

3, 3.5, 2, 5, 2 

3 

2 

4, 3, 4, 4 

4, 5, 3, 4 

4.5 - 5 

3 

2.5, 4, 7 

3, 3, 7 

4 - 4.5 

4 

3, 3.5, 4 

2.5, 3, 4 

6 

5 

4.5, 3.5, 4 

3.5, 2.5, 3 

7 

6 

5, 5, 4, 5 

5, 5, 4, 6 

5.5 

7 

6.5, 5, 7 

9.5, 9, 8 

8 

8 

2.5, 4, 4 

3.5, 3, 3 

2.5 - 3 

9 

3.5, 4, 3 

3, 3, 4 

3.5 

10 

3, 3, 3, 3, 4 

7, 3, 4, 4, 4 

6.5 

1! 

1, 3, 7 

2, 2, 7 

2 

12 

2.5, 5, 4 

7, 3.5, 4 

4 - 4.5 

13 

6, 7, 7 

8, 7, 6 

6.5 

19 

6.2, 5, 7 

8, 6, 7 

8 

20 

7, 4.5, 6 

7, 6, 5 

8 

14 

7, 3, 7 

8, 8, 8 

5 

15 

3, 3, 7 

6, 4, 7 

6 

16 

2.5, 4, 4 

3, 3, 3 

3 

17 

10, 9, 9 

10, 10, 10 

6 

18 

6, 8, 9 

8, 10, 10 

6.5 

21 

3, 3, 5 

2.5, 3, 3 

3.5 

22 

9, 9, 7 

10, 10, 10 

9 

23 

10, 9, 7 

10, 10, 7 

6 

24 

9, 5, 8 

10, 9, 9 

7 

25 

7, 4, 6 

10, 3, 6 

8 

26 

10, 10, 8, 8, 10 

10, 10, 10, 9, 10 

9 

27 

6, 5, 5 

9, 8, 7 '■ 

5 
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predicted value which should be expected for each test case. These 
values have been determined by interpolation between the pilot opinion 
isocontours of Figures 8, 9, and 10 which locate the test cases in terms 
of their short-period and Dutch Roll handling qualities, The agreement 
between ratings of overall handling qualities made by the pilot with 
the predicted values of each case is good for some cases but not for 
others. The primary factor of concern here is simulator realism. The 
tests were run with no visual cue given to the pilot, and a limited 
amount of kinesthetic cues, particularly in the normal direction, due 
to the simulator design limits and the washout system. The average of 
the overall pilot ratings for the three test runs are plotted along with 
the established contours for short-period and Dutch Roll handling 
qualities in Figures 27, 28, and 29. The dotted lines on these figures 
Indicate the handling qualities boundaries as suggested by the limited 
data taken in this study. While there is not exact agreement between 
the existing boundaries and the test data, the general trends exhibited 
by the two are in good agreement. 

For the short-period mode, the cases with good handling qualities 
(ratings of 4 and less) were not as close to the predicted values as 
were the cases where the handling qualities were not good (ratings of 
5 or more). For the Dutch Roll mode, the recorded ratings agree well 
with the predicted ratings for low Cooper-Harper ratings (below 4) and 
high Cooper-Harper ratings (above 7); however, agreement for marginal 
ratings (between 4 and 7) was poor. 
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Also, there appeared to be a qualitative relationship between 
control activity and pilot opinion of handling qualities. Figure 30 
depicts the variation of elevator activity for typical cases in which 
the longitudinal short-period characteristics were varied. As is 
evidenced from these traces, average Cooper-Harper rating of handling 
qualities increased as greater elevator activity was required to 
maintain level flight. The activity of the other control surfaces 
remained essentially constant while the short-period characteristics 
were varied. Similarly, variation of the Dutch Roll mode character- 
istics (Figure 31) were accompanied by changes in aileron activity. 
Higher Cooper-Harper ratings were accompanied by increased aileron 
activity necessary to maintain straight and level flight. The 
apparent activity levels of the remaining control surfaces were 
constant as Dutch Roll characteristics were varied. 

The only quantitative measures of control activity obtained were 
rms values of the control surface deflections. This measure was 
found to be unsatisfactory for determining control activity, due to 
the time averaging process involved in computing an rms of a parameter 
which is very small or zero for a substantial portion of the averaging 
period. 

Also, in the case of elevator rms deflection, it was not possible 
to exactly determine the equilibrium elevator deflection and bias this 
value out of the rms calculation. Another measure of control activity 
should be used in the future to quantify the activity level. Such 
quantities as the percentage of time the control surface deflection 
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exceeds a specified value or the number of peaks (max or min) per unit 
time experienced by a control surface may be more meaningful than rms 
for defining control activity. 

Taking the latter measure and u65ng a 1.5° threshold for both the 
elevator and ailerons, the total number of exceedances was determined. 
Figure .32 shows how elevator activity varied as the short-period 
handling qualities varied. As seen, this measure of pilot activity 
increases with undamped natural frequency and is not complimentary 
with the handling qualities contours. The t rend implies that the pilot 
uses more elevator motion to maintain straight and level flight when 
the short-period frequencies were high, possibly because he is not 
physiologically able to react rapidly enough to control a disturbance. 

In a similar manner, the number of times the aileron exceeded 1.5° in 
magnitude is plotted in Figures 33 and 3^ versus the Dutch Roll handling 
qualities. In this case, this measure of pilot activity is directly 
related to handling qualities and the region are complimentary; that is, 
pilot opinion of a configuration is worse when more aileron activity is 
required. 
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CHAPTER V 


SUMMARY AND RECOMMENDATIONS 

The interaction between ride and handling qualities for a typical 
STOL aircraft has been investigated. Simultaneous measurements of 
ride and handling qualities were made using a motion-based aircraft 
simulator and varying the short-period and Dutch Roll characteristics 
of the basic design. After plotting the ride quality associated with 
each test configuration along with the present requirements for short- 
period and Dutch Roll handling qualities, several important trends 
have been found. There appears to be a distinct trade-off between 
longitudinal short-period handling qualities and ride quality. Good 
Dutch Roll handling qualities and ride quality seem to be somewhat 
complimentary, but are, nevertheless, strongly related. It is up to 
the designers of future aircraft, particularly low-wing loading 
aircraft, to be aware of a vehicle's ride qualities as well as its 
handling qualities, and the results presented here should begin to 
guide him in these areas. 

Future work is recommended to investigate these interrelation- 
ships in a much more comprehensive manner. More test cases should 
be chosen to obtain measurements to better define how ride qualities 
vary as handling qualities are varied over their acceptable regions. 
Other stability derivatives than the ones used in this study could 
be chosen to vary the short-period and Dutch Roll handling character- 
istics to determine if individual derivatives have unique effects on 


73 


the levels of ride quality measured. This is particularly important 
to investigate for the lateral Dutch Roll mode. 

Likewise, the effects of variations in wing loading, turbulence 
level and power spectra, and flight conditions should be determined 
by repeating the tests as these additional parameters are changed. 
Ultimately, a general relationship could be developed to define the 
levels of ride quality to be expected of a given configuration as a 
function of all the above parameters. 


APPENDIX A 


UNIVERSITY OF VIRGINIA ANALOG FLIGHT SIMULATOR 

The simulator used in the preliminary studies consists of a fixed- 
base cockpit with rudder, stick, and throttle controls. The cockpit 
instrumentation includes artificial horizon, altimeter, rate of climb 
indicator, airspeed indicator, heading indicator (compass), and glide 
slope and localizer indications (ILS). The cockpit is connected to two 
analog computers, an EAI-580 and Pace TR-20, programmed with the six- 
degree-of -freedom equations of motion, ILS equations, and equations to 
compute the sum^of-squares of the lateral and vertical accelerations. 

A photograph of the instrumentation and the entire apparatus is 
presented as Figure 1 in the text. A more thorough description of 
the equipment may be found in references (25) and ( 26 ). 


APPENDIX B 


DERIVATION OF ANALOG EQUATIONS OF MOTION 

The s ix-degree -of -freedom equations of motion in Eulerian axes for 
a rigid body with an x-z plane of mirror symmetry are as follows (27): 

U + Q,W-RV + gsin©=2 X/m = a 

V + RU-PW-gsin$> cos 0 = E Y/m = a^ 

W + PV - QU - g cos $ cos 0 = E Z/m = a 


I l-l yl 

P - (R + PQ.) ~ + RQ (- ? - — ■£) = T- 


x 


I 


I 


X 


a + (P 2 . R Z) ^ w (LyJi) « m 

y y y 

R + (OR - P) ~ + QP ( X~* ) = 


I. 


I. 


Using the standard perturbation technique (24), the following 
nondimens ional equations are obtained: 


C Sq'u + C Sq a + ~ C Sq 0 + — C Sq a + C Sq 6 
x 0 ° x 0 ° 2U n x oo 2U n x. oo^ x co e 
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Aa I u 
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These are rearranged for ease in analog programming. An order of 
magnitude analysis allows the products of inertia and differences in 
moments of inertia terms to be neglected compared to other terms in 
the moment equations, yielding 
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APPENDIX D 


SIGN CONVENTIONS 

Sign conventions for forces, moments, angular rates, and linear 
velocities are shown in the sketch below. The positive sense of each 
variable is shown. 



Sign conventions for control deflections are listed below: 

Control Surface Positive Deflection 

Elevator Trailing edge down 

Ailerons Left aileron trailing edge down 

Right aileron trailing edge up 

Trailing edge left 


Rudder 



APPENDIX E 


DERIVATION OF EXPRESSION FOR NORMAL ACCELERATION 
DUE TO TURBULENCE 

If one eliminates the force equation in the X body axis direction, 
assumes variations in u to be negligible, neglects C 2 and C z and 

q & 

assumes, a sinusoidal forcing function, the following equations 
result (24): 


w 

(2ys - C z )a - (2y)q = (-C z ) 

a a 0 


-< C *. 5 + C r> + (, B 5 - C m “ <' C 

a a q 


a 


w 

C s + sC ) 7 . 2 
m. m ' U ft 

a q 0 


where y - m/pS£ 

» B c yp s * : 


SL - c/2 

s = Laplace operator 
w = gust input magnitude, 


One may solve these expressions for ot/w and q/w Q by Cramer's rule 

37 37 


to find 


■s— = ■ "■t [ ( — C “ C s + C s) (2y) ■ C ( i pS - C )] 

w_ U«A m^ m, ■ ■ z_. ' B m. 


9 0 


a 


3- = — [ (2Us - C ) (-C - C s + C s) - C (C sC ) ] 

w U,A z ' ' m m, m„ z m. m 

lx a a q a a a 


PRECEDING page blank NO? mMm 


A - (2ps -C )0 B . - C ) -2u(C m _s + C m ) . 

a q a a 


The transfer function for the acceleration factor is given by 




w gc w w 

g g g 


Making the substitutions for q/w g and a/w^, the following expression 
results: 


5- = -=% (-2C C ) s . 
W g 9 cA m , Z a 


Setting s = ioj, one may form the expression: 


2 kU* 
n i 0 


4C 2 C 2 w 2 
m z 
q a 


W g g 2 c 2 [ (“2yC -1 R C -2yc )V+(C C -2pC -2y i R a> 2 ) 2 ] 

m q b z a z a m q m a B 

As a measure of the total normal acceleration, we wjl 1 form 


-r is-i • 


‘N > <w 

_oo g 


This may now be written as: 


2 

fOO U) ckl) 

*N ~ K j 2 2 

—ob A + Boj +oj 


C C - 2uC 

where z m H m 2 

A-( ' %Mj B g ) 8 " d 


- *M*9*»«^ 




(-2WC - i t - 2UC r - i (C C - 2|iC ) 

m q 6 2 o m d 8 2 c m q m a 


Following the development In Reference (24) for an approximation 

of the short-period mode, the following relations are obtained for 

a) and S 
n s 5 

c C - 2yc 

.. 2 2 a m q m a 


2PC + i_C + 2yC 
m ■ B zl, r 


2C w 
s n 


By use of the two previous expressions for m and 2c u) , it can be 

s s n s 

seen that 


A = a) 


Solving 


B = 20) n 2 (2? s 2 - 1) . 

s 


oi 2 dw 


A + Bw + to 


we attempt to write this in the form 


2 

a) dw 


(E + w 2 ) (F + to 2 ) 


where 


E.F = 


B ± YB 2 - 4A 



We now examine the quantity B - 4A: 


B 2 '- •* - [2o> 2 (2C s 2 - ))] 2 * 4 - S, 4 C S 2 <C S 2 - ') • 

„ s s s 


From this, we may identify three different cases for the roots of the 
characteristic denominator. 

Case 1 . C s ^ < 1 


Complex roots exist since - 4A < 0. Integrating 


r 


2 

0) dw 


2 5 

—co a + Bto + io 


over the complex plane, we obtain the following result after much 
complex arithmetic 


2 . 

a) dw 


r . v— 

A + Boo 2 + 0 )^ 2(0 “ S 2 


Case 2. C s = 1 


Two real equal roots exist since B z - 4A =0. Integrating, we find 


2 

!°° (0 d03 _ 7T 

/ 2 ~ 2u) 

- 00 A + Boo + 0) n 


Case 3. C s ” > 1 


Two real unequal roots exist since - 4A > 0, Again, integrating 
we find 




J°° co da) _ TrjE/F - Ft/E) 

— » A + Baj^ +w (E - F) 

»|2; s 2 - 1 + z; s V? 5 2 -nt2g s 2 - I - z^VeT-il* 5 

S 

[2C s 2 - 1 - *-]][&* - 1 + 

C a/s -1 

n_ s ? s 


If one plots the numerator function over the range 1.1 < C g < 2.0, 
one finds It nearly linear in ? ■. We may then write 


w 2 dco 


J 2 U 

- 00 A + Boj +o) 




In the original expression for a N> a constant K appeared which was 

2 2 

a function of the product C . C .In this study C is held fixed 

z m z 

o q a 

and C is used as a variable to modify the short-period character- 

m q 

istics. An examination of test cases 1, 2, and 13 (C = 1.9098) 

a 

reveals that there is almost a linear relationship between C and £ . 

m q 5 

A mathematical fit to this trend is r = -0.011 C + 0.332 or 

s m 

q 

C ~ 33(1 - 3? s ). Incorporating this function into the three 

q 

expressions gives 




w 

If 
1 ■ 


\y 

H' 

f 

lV 

P' 

U: 






5 S < 1 

5s-’ 

5s > ’• 

where the constants Kj , K^, anc! are functions of flight conditions, 
aircraft mass, pitch moment of inertia, and the stability derivative 

V 

Expanding the radicals by the binomial expansions yields 

U/V^7) = 1 - 2 ?2 + 8 ^ + H -°' T * 

0/V? 2 - 0 = + | J 5 + h.o.t. 

c c 

Thus 

0 - 3; ) 2 , 

■ - 1 - 6c; s + 8.5;/ + H.O.T. ; s < 1 

V 1 - 5s 2 

(I - 3; ) 2 , 

— — - 9C S - 6 - £ * H.O.T. 5, > 1 

5s - ' 



a N * K 1 


"gC - 3 5 ,)■ 

v/ 1 5s' 
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a = K_ —2— 

N 2 a) 

n s 


a - K 
a N K 3 


“n Vc - 1 
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and simplified expressions for a^ are 


N 


1 0 ) 


n 


w 

- K -2— 
2 “n 

s 





- 6 - - 3 - + H.O.T.) 

2 5s 



it should be noted here that the mathematical fit for C versus 

m q 

r Is unique to the particular value of C chosen. However, fits for 
s m 

a 

other values of C = constant exhibited the same trends as above, 
m 

a 

that is C ~ (1 - kr ) where k is a positive constant. Therefore, 
m s 

q 

in general 
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